Cabibbo allowed D Kn^ decays 

S. Fajfer"^^, A. Prapotnik^ and P. Singer*^ 

a) Department of Physics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, 

Slovenia 

b) J. Stefan Institute, Jamova 39, P. O. Box 300, 1001 Ljubljana, Slovenia 
c) Department of Physics, Technion - Israel Institute of Technology, Haifa 32000, Israel 



ABSTRACT 

The weak radiative Cabibbo allowed decays —>■ K^n^'j and —>■ K^tt^'j with 
nonresonant Kn are investigated by relying on the factorization approximation for the 
nonleptonic weak transitions and the model which combines the heavy quark effective 
theory and the chiral Lagrangian approach. The dominant contributions to the amplitudes 
come from the long distance effects. The decay amplitude has both parity violating 
and parity conserving parts. The parity violating part includes also a bremsstrahlung 
contribution. The branching ratio obtained for the parity conserving part is of the order 
10~^ for the — > K~7i~^'-f decay and 10~^ for —>■ K^n'^'j, when the effect of light 
vector mesons is included, and smaller otherwise. The branching ratio for the parity 
violating part with a photon energy cut of 50 MeV, is close to 10~^ for the D° decay 
and 4 x 10~^ for the decay. We present Dalitz plots and energy spectra for both 
transitions as derived from our model and we probe the role of the light vector mesons in 
these decays. 



I. INTRODUCTION 



The investigation of radiative and dilepton weak decays of pseudoscalar charm mesons 
has been pursued rather vigorously in recent years, both theoretically and experimentally. 
To a certain extent, this activity has been fueled by the ongoing search for physics beyond 
the standard model, which might be of measurable consequence in certain charm radiative 
and dilepton decays [1-4]. To date, no radiative or dilepton weak decay of D has been 
detected. However, upper bounds have been established for a sizable number of these 
decays. The radiative decays p°, uo^ ^ 0, K*^ +7 were recently bounded P| to 

branching ratios in the 10~^ range, which is approaching the standard model expectations 
(see, e.g. [^, 0] where additional previous works are mentioned). The dilepton decays 
D — > Pl~^l~ , D — > Vl~^l~ are the subject of intensive searches at CLEG and Fermilab 
Here again, with upper bounds of 10~^ — 10~^ for branching ratios of the various modes 
one approaches the expectations of the standard model [|, |], H, The situation should 
improve in the future, due to new possibilities for observation of charm meson decays at 
BELLE, BABAR and Tevatron. Recently, upper limits in the 10~^ — 10~^ range were 
established |T0| also for D° dilepton decays with two nonresonant pseudoscalar mesons 



in the final state D'^ —>■ {tt^7t~ , K'n'^ , K~^K~) though no comparable results are 

available yet for similar photonic decays. 

In the present work, we undertake the study of the Cabibbo allowed radiative decays 
—>■ K^n^'-f and K^n'^'j, which we consider to be the most likely candidates 

for early detection. These decays are the charm sector counterpart of the K mi'y 
|TT|, |r^, |T3] decays, which have provided a wealth of information on meson dynamics. 
In the strange sector, the 7r+7r°7 and — > 7r+7r^7 are singled out as the most 

suitable ones for the investigation of the radiative decay mechanism; this, since the relative 
suppresion of the corresponding — >■ vr+vr^ and — ^ tt'^tt' amplitudes leads to a 
situation where the direct radiative transition is not overwhelmed by the bremsstrahlung 
part. 



In the K nnj decays, the long - distance contribution is dominant |13]. In the 



charm radiative decays, the theoretical studies show that likewise, the long - distance 
is the dominant feature of the decays 0, ^ ^ [7], ^. The short - distance contribution 
realized by the penguin diagram c ^ wy fl^, |15], [1^ might play a role in certain Cabibbo 
suppressed decays, which are not discussed in the present paper. 

In the charm sector, the nonleptonic D meson decays still provide a continuing theo- 
retical challenge (see, e.g. p!7| , |T8| and references therein). The short distance effects are 
considered well understood but the perturbative techniques required for the evaluation 
of certain matrix elements are based on approximate models. Usually the factorization 
approximation is used (see, e.g. [|I7|, |I^), although the experimental data indicate the 
apparent need for the inclusion of nonfactorizable amplitudes in certain channels. 

In this first treatment of -D — Kn'y decays we use the factorization approximation 
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for the calculation of weak transition elements. We consider the use of this approach to 
be justified by the "near" success of the approach for the nonleptonic amplitudes. This 
will involve its use in the D{D*)Ktt vertices as well as in D{D*) V and D{D*) —>■ P 
transitions, all of which required for the calculations of the D -—>■ Kn'-f amplitudes within 
our model. For the evaluation of the {D,D*) — >■ (P, V) transitions, we use the information 
obtained for these matrix elements from semileptonic decays (see, e.g. The general 

theoretical framework for our calculation is that of the heavy quark chiral Lagrangian 



21| , p2| . In the K 7r7r7 decays, it has been shown that intermediate light vector 
mesons play an important role in the decay amplitude We shall investigate the 

role of intermediate light vector mesons also in the D Ktt'j amplitude. In order to 
accomplish this, we use the extension of the formalism of [^, ^ to include also the light 
vector mesons p3|, p4| . 

The present study of K^n^'-f and K^n^'-f shows that the direct part of 

the radiative amplitude is not much smaller in strength than the bremsstrahlung part, 
rather similarly to the inhibited K decays mentioned above. If confirmed by experiments, 
it places these decays in the status of a most suitable ground for the investigation of the 
mechanisms involved in such nonleptonic D decays. 

In Section II we present the theoretical framework for our calculation. In Section III 
we display the explicit expressions of all the calculated decay amplitudes. Section IV 
contains the discussion and the summary. 



II. THE THEORETICAL FRAMEWORK 

The nonradiative two-body D - decays, from which the bremsstrahlung part of the 
radiative decays originates, are K^ir^ and — > K^Tr~^. The weak AI = 1 

transition leads to two independent isospin amplitudes in the final state, A1/2 and A3/2 
and the relations to the physical decays is p 



A{D+ ^ K'n+) = A3/2; K-TI+) = + ^^3/2. (1) 

From the determined branching ratios |^ of BR{D^ K^n^) = (2.89 ± 0.26)% and 



BR{D^ —>■ K~7T^) = (3.83 ±0.09)% one learns that the relative size of the absolute values 
of the amplitudes \ A(D'^ — > K~ 7r~^)\/\ A(D^ K^7r~^) \ is 1.84. Using also the information 
from the third decay, A{D'^ K'^n^) = ^^1/2 — -^^3/2) the isospin analysis shows that 
1^1/21/1^3/2! — 2.7, and their relative phase is 90° p6|]. Despite this knowledge, there is 



still no complete interpretation for the mechanisms leading to the decays although 
it is clear that the situation is different from the K Tin channels, where AI = 1/2 
enhancement introduces a large disparity between the final state isospin amplitudes. The 
relevance of the above picture to the radiative decays will be discussed in the last Section. 
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Since our problem of describing the D Kir^y decays involves transitions between 



heavy mesons and light pseudoscalars, we adopt the effective Lagrangian |21|, ^ which 
contains both the heavy flavor and the SU {3)l x SU{3)ii chiral symmetry as the theoretical 
framework for our calculation. From the experience with K — > mr'j decays, one knows 



TT| , |T2|, |T3[ that the decay amplitude is largely determined by contributions from virtual 
vector mesons. Considering the possibilty that vector mesons would play a role in the 
D — s> KiT'y decays as well (we remind the reader that we consider here nonresonant Kn, 
the decays D K*'y having been treated separately |16|, we should complement the 



Lagrangian by introducing light vector mesons. For this we choose the generalization of 
the original Lagrangian ^ by Casalbuoni et al ||2^ in which the original symmetry 
is broken spontaneously to diagonal SU{3)v with the introduction of the light vector 



mesons. We present here in some detail this formalism (for more details see ^^), which 
we use as the main tool of our calculation. We shall perform the calculations also without 
vector mesons in the Lagrangian in which case the original heavy quark chiral Lagrangian 



21| , |22[| is used, in order to clarify their role in these decays. 

The light degrees of freedom are described by the 3x3 Hermitian matrices 



n 



, VS I VP 

7l~ 



V 



and 



K- 



VS 



VO 



-7r" I Vs_ J 



K+ 



(2) 



V2 

Pm 



p 



V2 



(3) 



for the pseudoscalar and vector mesons, respectively. They are usually expressed through 
the combinations 



u = exp f ^ ) , (4) 



/ 

where f '2=l = 132 MeV is the pion pseudoscalar decay constant, and 



P/^ 



(5) 



where = 5.9 was fixed in the case of exact flavor symmetry In the following we 
will also use the gauge field tensor F^y{p) 
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Ff^uip) = d^p^ - dupf, + [pf,, pu]- (6) 

It is convenient to introduce two currents V^j = \{u^D^u + uD^u^) and = \{u^ D^^i^u — 
uD^u^). The covariant derivative of u and is defined as D^u = {d^ + B^)u and 
Dfj,u^ = {d^ + B^)u\ with = ieB^Q, Q = diag{2/3, —1/3, —1/3), 5^ being the photon 
field. 



The hght meson part of the strong Lagrangian can be written as [ p8[] 
Cught = -y{tr(^^^^) + atr[(V^-p^)2]} 

The constant a in (|^) is in principle a free parameter. In the case of exact vector meson 
dominance (VDM) a = 2 |2^. However, the photo production and decays data 
indicate that the SU (3) breaking modifies the VDM in 

U-y = -egj'B.ip'^ + - ^<D'^). (8) 
Instead of the exact SU{3) limit {g^ = mv/ f), we shall use the measured values, defining 

<Viey,q)\V,\0>=te;iq)gyiq'). (9) 

The couplings gv{my) are obtained from the leptonic decays of these mesons. In our 
calculation we use 5'p(m^) ~ gp{0) = O.nGeV"^, g^jim"^) — g^i^) = O.lbGeV"^ and 
g^iml) g^{0) = 0.2AGeV\ 

Both the heavy pseudoscalar and the heavy vector mesons are incorporated in a 4 x 4 
matrix 

Ha = 1{1 + MP:,Y - Pal,) , (10) 

where a = 1,2,3 is the SU{3)v index of the light fiavors, and P*^, Pa, annihilate a spin 
1 and spin heavy meson Qqa of velocity v, respectively. They have a mass dimension 
3/2 instead of the usual 1, so that the Lagrangian is in the heavy quark limit ttiq — > oo 
explicitly mass independent. Defining moreover 

Ha = i'Hy = {p:iY + pli,)l{i+i^), (11) 



we can write the strong Lagrangian as 



Ceven = Cught + iTr{HaV ^D^ Ha) 

+ igTr[H,-f^-f,{A^),aHa] + iPTr[H,v^{y^ - p,)taHa] , (12) 
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where D^'Ha = {d^ + - ieQ'B^)Ha, with Q' = 2/3 for c quark. 

The couphng g can be fixed 0] by using the data on D* Dtt decay width. 
These data give g = 0.59. The plus sign is taken to be in agreement with the quark model 
studies. The parameter P is less known, but it seems that it can be safely neglected 

The electromagnetic field can couple to the mesons also through the anomalous in- 
teraction; i.e., through the odd parity Lagrangian. The contributions to this Lagrangian 
arise from terms of the Wess - Zumino - Witten kind, given by pOl 



4i = -4^6^-^Tr(9,p.9„p^n). (13) 
The coupling Cyvn can be determined in the case of the exact SU{'i) fiavor symme- 



try following the hidden symmetry approach of j^S], ^] and it is found to be Cyvu = 
3gl/32TT'^ = 0.33. In the actual calculation, we allowed for SU{3) symmetry breaking and 
we used the VP'j coupling as determined from experiment p5|. We will also need the 



odd - parity Lagrangian in the heavy sector. Such terms are required by the D* — > D'y 
transition, which cannot be generated from (0). There are two contributions P^ , [ 
it, characterized by coupling strengths A and A'. The first is given by 



m 



£i = t\Tr[Haa^,F^''ip)abH,]. (14) 

In this term the interactions of light vector mesons with heavy pseudoscalar or heavy 
vector mesons is described. The light vector meson can then couple to the photon by 
the standard VDM prescription. This term is of the order 1/ A^ with A-^ being the chiral 



perturbation theory scale [35 



The second term gives the direct heavy quark - photon interaction and is generated 
by the Lagrangian 

£2 = ->!Tr[Haa,,F^\B)Ha]. (15) 

The parameter A' is given in heavy quark symmetry limit by A' ~ —l/{Qm^ and it 
should be considered as a higher order term in 1/ itlq expansion . 

In order to gain information on these couplings one has to use the existing data 
on D*^ — > -D°7, -D*"*" D^'-y and D^^ — > D^'-y decays. Experimentally, the ratios 
i?o = r(D*o ^ L)°7)/r(D*o ^ Z^o^O) and R+ = T{D*+ D+-f)/T{D*+ D+n^) are 
known ||2^. These data determine two possibilities |2^. One of them is \\/g\ = 0.839 
GeV~^, \y/g\ = 0.175 GeV~^. The second one does not agree with present data. With 
g = 0.59 we obtain A = ±0.49 GeV'^ and A' = ±0.102 GeV-\ 

The A' ~ —l/lQuic) would give with the mass of charm quark rric = 1.4 GeV that 
A' = —0.12 GeV^^, in good agreement with the above value. The simple quark model 



analysis indicates that A' and A are both negative |]36[. In our numerical calculations we 
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give the results using these parameters. In the hterature (e.g. instead of 

A the (3 parameter is often used. The value (3 = 2.3 GeV~^ corresponds to A = —0.49 
GeV-\ since 2A^ {9p/ml + gj{3ml)) = -{2/3)(3. 

In addition to strong and electromagnetic interactions, we have to specify the weak 
one. The nonleptonic weak Lagrangian on the quark level for the Cabibbo allowed decays 



can be written as usual 19 



CjiliAc = As = 1) - 

where Oi = {ud)Y_^{sc)v-A,fi and O2 
elements, Gp is the Fermi constant and {^i'^^ 



G p 

'71 



KdV;:[aiOi + a202] 



(16) 



{uc)v-A,t,{sd)'^_j^, 



Vij are the CKM matrix 
i'^ = 'I'i7''(l — 7^)\l/2. In our calculation 
we use ai = 1.26 and 02 = —0.55 as found in |T9[| . 

At the hadronic level, the weak current transforms as (3/,, 1r) under chiral SU{3)l x 
SU{3)ii, is linear in the heavy meson fields and D*" and is taken as 



Jo": 



\ iaTr[Y{l-^,)H,ul] 

+ «2Tr[7^75i/bt;«(p"-V")fc,My + ... 



(17) 



where a = fny/^H [^|, «i was first introduced by Casalbuoni et al. while 0:2 was 
introduced in [|^. It has to be included, since it is of the same order in the l/mg and 
chiral expansion as the term proportional to ai |^0[ . 

The relevant matrix element is parametrized usually in D Vlvi semileptonic decay as 



16, 19, 20, 38 



rriD + my 



2mv 



+^^v-q—^qtJi{M{q ) - ^o(g )) + i{mD + mv)[tv^Ai{q 



i{p + Pv)^.A2{q^)] 



;i8) 



mp, + mv 

where q = p — pv- In order that these matrix elements should be finite at = 0, the 



form factors satisfy the relation |]T9 



TUh + my 



2my 2mv 
and ^3(0) = Aq{0). We take the following expressions for the form factors at g^^^ 



(19) 



M + m 
'M + A'' 



(20) 



2 ^ —V^aiQ^y/M 



and 



MiL.) = ■ (21) 

2 N 2g^M + m 

Mimax) = ^^^^ 

where A stands for the and D mass difference. Assuming the pole dominance one 
can connect the value of form factors at q'^^x momentum transfer by F{0) = 
F{qmax) (l - (M - my/M^y where F stands for V, Ai or M is the D meson mass 
and m is the light vector meson mass. Using the experimental data l^'^'^*(0)l = 
1.02±0.12, |Af^*(0)| = 0.55±0.03 and |Af^*(0)| = 0.40±0.07, we find for the couplings 
A = -0.56 GeV-\ \ai\ = 0.156 GeV^/^ jasl = 0.052GeV^/l The value of A is in good 
agreement with results obtained from D* D'y data. 
The light weak current is derived to be 



Jt, = tf{u[A^ - a(V^' - p^]u%. (23) 

The photon emission is obtained by gauging the weak sector too. The important con- 
sequence of this procedure is that thereby the gauge invariance of the whole amphtude 
is achieved. This turns out to be equivalent to the usual procedure of achieving gauge 
invariance in bremsstrahlung processes with a momentum dependent strong vertex, as 
pointed out for the somewhat similar process V — > PP''y. Actually by gauging the 



weak sector we produce the same graphs, which were necessary to induce to satisfy the 
gauge invariance pQ . 



III. THE DECAY AMPLITUDES 



The general Lorentz decomposition of the D{P) K{p)iT{q)'-f{k,e) decay amplitude 
is given by 

+ F^e^'^^^e^v^kf.q,^ . (24) 

The part of the amplitude containing the Fi form factor is parity violating, while the one 
with F2 is parity conserving. Both of them are functions of scalar products of momenta 
a.s k ■ p, k ■ q. Note that Fi contains contributions which arise from bremsstrahlung part of 
the amplitude as well as a direct electric transition. On the other hand, F2 corresponds 
to the magnetic transition. 

In order to determine Fi, F2 we use the model described in the previous Section. The 
diagrams contributing to these form factors are given in Figures 1-4. In Figures 1 and 
2 are given Feynman diagrams contributing to the K^tt'^'j decay amplitude while 

the contributions to the K^-k^'-j decay amplitude are presented in Figures 3 and 4. 

Note that we denote heavy mesons by one full line, light pseudoscalar mesons by dashed 



q ■ e p ■ e 
q ■ k p ■ k 
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lines, light vector mesons by two full lines and photons by wavy lines. The weak vertex 
is denoted by square box. 

Before proceeding to the actual calculation, we note the following complication. As 
well known, the leading terms of the expansion of the radiative amplitude in the photon 
momentum (k) are determined [^^j by the original amplitude {D — >■ Kir in our case). 



However, the nonleptonic D Ktt amplitude cannot be calculated accurately in the 
factorization approximation from the diagrams provided by our model. Such a calculation 
gives a rather good result for the K^n^ channel but is less successful for the 

K'^TT~ decay. In order to overcome this defficiency and to be able to present accurately 
the bremsstrahlung component of the radiative transition, we shall use an alternative 
approach for its derivation. This approach then is to use the values of the experimental 
amplitudes D —>■ Ktt, assumed to have no internal structure, for the calculation of the 
bremsstrahlung component. In order to accomodate this we rewrite the decay amplitude 
(24) as 



q ■ e p ■ e 
q ■ k p ■ k 



+ F, [{q .e){p.k)-{p- e){q ■ k)] + F^e^''^^ e^v^k^q^ 

(25)' 

where Fq is the experimentally determined D Kit amplitude and Fi, F2 are the form 
factors of the electric and magnetic direct transitions which we calculate with our model. 
When intermediate states appear to be on the mass shell, we use Breit Wigner formula. 
However, we remark at this point that since we are interested in the D —>■ {KTT)nonres1 
transitions, we delete the region of the K* resonance appearing in diagram D^2 
retain only the region in {p + q)"^ which is beyond rriK* ± Tk*- 

In Appendix A we present explicitly expressions for the form factors for the decay 
K^n^'j using the notations Af, Bf, etc. . The amplitude Af , Bf, etc. is 
obtained as a sum of the amplitudes presented by the corresponding diagrams in the i-ih 
row in Figures 1-4. Each Af (or Bf, etc) is gauge invariant. For the electric parity - 
violating transition, we define both the total amplitude provided by the model Fi, as well 
as the direct part only, Fi, obtained after deleting the bremsstrahlung diagrams. Then 
the amplitudes for K^n'^'y are 

Fi(D+ K\+^) = jz^At + C+), (26) 

1=1 

F,(D^^KV-,)^^-±-^^j:(M^Ct), (27) 

F^(D* /?%+7) = Y.[Br + Dt). (28) 

1=1 



8 



In the case of ^ K 7r+7 we have 

F,{D'^K-7r+^)=j:{A'^ + C% (29) 



i=l 



F,(D'^K-.^,)^^-±-^^Y.Al (30) 
F^(D'' K-7,*',) = + (31) 

i=l 

where etc. are gauge invariant sums of the amphtudes arising from the graphs in 

the 2-th row. In Appendix B we present the form factors for the —>■ K'tt^'j decay. 
We denoted by B^'^ contributions which are created by Oi operator and by C^'^, 

D^'^ contributions caused by 02- 

As we mentioned, in the calculation we used the experimental value of A{D Kit) to 
calculate the bremsstrahlung part; Fi is then calculated by subtracting the bremsstrahlung 
component from Fi. 

The differential cross section of the decays is given by 

= p^^i-^i'*"-*"-- P2) 

where Ai is the decay amplitude, given by Eqs. (|^) or (^), = (-f ~ ^)^) 
"^23 — ~ pY-i where P, p are respective four - momenta of D- meson, photon and 
K meson. The total decay width is a sum of the parity conserving expressed by F2 and 
the parity violating contributions expressed by Fq, Fi, T = Tpc + Tpv (the PC and PV 
amplitudes do not interfere in the total width). Before giving numerical results, we make 
a few comments. 

The expressions for the amplitudes given in the Appendices contain several constants. 
A few are well determined (we use values given in [^) and require no further explanation; 



as to the rest, for Jd we use the lattice result, Jd = 207 MeV |^ and for f^* = l.lSfo- 
The couplings g, A, A' are determined as previously explained and we use g = 0.59, 
A = -0.49 GeV"^ and A' = -0.102 GeV-^ 

Some of the amplitudes, like A21, ^245 ^255 ^tc, contain the weak transition 

D* to TT. This requires the vr to have very large momentum, which means that such 
graphs, which vanish in the heavy quark limit, give very small contributions. 

Turning now to the presentation of the results we have to start with a discussion of 
the bremsstrahlung contribution (IB). In our model IB is given by diagrams (IB)" = 
i + ^2 i + ^1 i) for the -D° — > K^n^'y decay (the first two rows and the fourth row of 
Fig. 3) and by diagrams (IB)^ = + ^tj + + C'2'^j), i-e. the first two rows and 

rows 5 and 6 of Fig. 1 for the — > K^n'^'-f decay. Now, in the limit of vanishing photon 
energy, the first two terms in the expansion of the IB amplitude in terms of the photon 
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energy, obey the Low theorem [^. Ahhough this is fulfilled theoretically, the question 
arises wheather the D —>■ Ktt amplitude, as derived from our model, describes correctly 
the observed —>■ K^n^, K~'7r'^ decays. We calculated the amplitudes of these 

decays using our model and we find that the branching ratios obtained with factorization 
approximation are 4.1% and 17% respectively, compared with observed branching ratios 
of 2.9% and 3.9% [^. It appears that although the model is reasonable for — »• K^ir'^ 
(the AI = 3/2 amplitude), it misses the amplitude of K^n^ by a factor of 2. 

On the one hand, this gives us a certain reassurance on the suitabilty of the model 
we use for calculating the radiative amplitudes. On the other hand, we shall perform 
also an alternative calculation, whereby the bremsstrahlung amplitudes of the model are 
deleted from the total radiative amplitude and replaced by the "experimental amplitude" . 
This procedure is undertaken in order to enforce the fulfilment of the Low theorem for 
our radiative amplitudes. Thus, we assume constant D Kn amplitude of correct 
magnitude to reproduce the observed rates of —>■ K^Tr~^ and —>■ K~tt~^, from which 
we calculate the bremsstrahlung (IB) amplitudes. These have the form of the first term 
in (^) with constant Fq. To this we add the F2 terms of the magnetic transition, which 
is not affected by this procedure, as well as the parity-violating Fi terms not belonging to 
(IB)^ and {IB)~^ diagrams. These Fi terms then represent the direct electric transition 
of the radiative amplitude. We present results for both these alternative procedures. 
Although, the procedure based on the experimental D Kn amplitudes is apparently 
more reliable, we consider the "model" calculation to be of intrinsic value, setting out the 
ground for future calculations. 

There is one more item to be explained. We are interested in the role played by the 
vector mesons in these decays; obviously not in the direct Kit channel, which belongs to 
_^*°7 and was treated separately [^, |^ |T6[, rather as they appear as intermediate 
particles in VDM (e.g. diagrams ^3^4, 6*3^4, -02^3, ^3^3 and others) or in the crossed channels 
(e.g. diagrams A^^i, D^i, D^2 C*?3) -^3^2; others). This is the main reason for our using 
an effective Lagrangian which contains the light vector mesons [EB] . However, we calculate 
the radiative transitions also without including vector mesons in the Lagrangian, i.e. we 
drop all diagrams containing a double line in Figs 1-4 (gauge invariance is maintained), 
which allows to elucidate their role in these decays. 

For the parity conserving part of the decays, representing the magnetic transition, we 
obtain 

BR{D+ K\+-f)pc = 2.0 X 10-^ (33) 

BR{D^ K-Tx+-i)pc = 1.4 X 10"^ (34) 

If we disregard the contribution of vector mesons, the rates are reduced to BR{D^ 
i^°7r+7)^°(V^^ = 3.0 X 10-6 and BR{D^ 7^-71+7)?,°^^ = 6.6 x 10"^. The decrease is 
sharper for the D*^ decay, since in this case the light vector mesons gave the dominant 
contribution to the rate, this is not the case for where such a contribution is doubly 
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Cabibbo suppressed. The differential distribution for these transitions, as a function of 
"^12 — ^ is given as the dashed hne distribution in Fig. 5a, 5b, for these two 
decays. The distribution is mainly symmetrical, with the peak occuring at A; ~ 400 MeV. 
Thus, this is the region in which the effect of the direct transition has best chance for 
detection. 

Turning to the parity - violating transitions we start with the procedure whereby we 
enforce the Low theorem by using Eq. (25). Here we face the question of unknown phase 
between Fq and Fi. We give therfore the results in terms of a range, limited by minimal 
and maximal interference between Fq and Fi. 

Thus, we get for the branching ratios of the electric transitions, with |Fo| determined 
experimentally, 

BR{D^ ^ i?°7r+7)^>?°f = (3.6 - 3.8) x 10-^ (35) 

BR{D+ K%+-i f>^^l^'^ = (2.3 - 2.5) x 10^^ (36) 
For the radiative decay we get 

BR{D^ ir-7r+7)|,>^°f = (5.0 - 15) x 10-^ (37) 

5i?(D° ^ ^"7r+7)|>|0°^^^ = (2.6 - 11) X 10-^ (38) 

The uncertainty in the Fq/Fi phase is less of a problem in K^n^'j than in — > 

K~7T^'y. If we take the bremsstrahlung amplitude alone as determined from the knowledge 
of \Fq\, disregarding the direct electric Fi term, the above numbers are replaced by 3.6 x 
10"^ and 2.3 x 10"^ for decay and 8.6 x 10"^ and 5.5 x 10"^ for the D'^ decay. 
In Fig. 6 we also show the dependence of the branching ratio of the bremsstrahlung 
amplitude on the lower energy bound, for both decays. The contribution of the direct 
parity violating term (putting Fq = ), is BR{D^ K'^n~^'y)iiir,pv = 1-0 x 10~^ and 
S/?(L)° ^ K-7r+^)dir,pv = 1-64 x 10"^ 

We checked also for the PV-transition the effect of the vector mesons. Using the 



Lagrangian of Refs. [^, ^ we found that in the PV-transition the effect of vector mesons 
is rather negligible; there is practically no change in (p5D, ( p6D and only a a narrowing of 
the range in (|37|) and (|38D , to bring it essentially to the values of pure bremsstrahlung we 
indicated after Eq. (^). 

We have calculated the decay rates also by using our model, for the whole radiative 
amplitudes i.e. using all graphs of Figs. 1-4. Comparing these results with those of Eq. 



]) - (p8D gives an indication of the possible uncertainty of our model. We obtain 

BRiD^ A^°vr+7)^>v?wer = 3.0 x 10-^ (39) 
BR{D^ ^ i^V^)^>yZlf = 2.5 x 10"^ (40) 
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For the radiative decay we get 

BRiD' K-7r^^)%%'::^ = 2.3 x 10-^ (41) 

BRiD' K-n^^)%Zlf = 1-5 x 10"^ (42) 

In Fig. 5c, 5d, we compare the rate of the decays dT = dTpc + dTpv for the two alter- 
native calculations concerning the PV part. In Fig. 5e, 5f we compare the rate of decay, 
dr = dTpc + dTpv calculated from (pSD to the bremsstrahlung rate, to emphasize the 
feasibilty of detecting the direct emission. Finally, in Fig. 7 we present Dalitz plots for 
these decays. 

IV. DISCUSSION AND SUMMARY 

The calculation we presented is the first attempt to formulate a theoretical framework 
for decays of type D KiT'y, with nonresonant Kit. The calculational framework is 
the strong Lagrangian (12)-(15) used in the tree approximation, and factorization for the 
weak matrix elements. In the present article we treat the Cabibbo allowed decays and 
among these, only channels which have both inner bremsstrahlung and direct radiation 
components, K^n^'-) and K~7r~^'y, i.e. those which are the most likely ones 

for early detection. There is a third channel in this class, — ^ K^tt^'j, which has only 
a direct component in the radiative decay and will be discussed separately. 

Our results show that the relative expected strengths of the direct and bremsstrahlung 
components are of a magnitude which would permit the experimental determination of 
both, in next generation of experiments. This is important, since the direct amplitude 
provides information on the decay mechanism. In the radiative decay of D~^, the magnetic 
direct component amounts to about 6% of the total rate (see Eqs. (33), (35)) and together 
with direct electric component which is of comparable magnitude (see Fig. 5a), dominate 
the decay spectrum in the region of high photon energies, say above k ~ 250 MeV (see 
Fig. 5f). A similar situation occurs in the K~'TT~^'~f transition, where the direct 

radiative decay containing both electric and magnetic parts which are of nearly equal 
magnitude, amounts to over 30% of the total radiative decay rate (see Eqs. (34) and 
(37)). The numbers we mentioned are for k > 50 MeV, but as we stressed the region of 
high photon momenta beyond 200 MeV is where the direct transition even dominates. 

The above large relative rates of direct/IB are somewhat similar to the occurences 
in K2 — >■ 7r+7r~7 and n^ir^'y, although in the D decays there is no apparant 

suppression of the original decays. This is most probably related to the mechanism of the 
decay provided by the heavy quark Lagrangian and the coupling constants involved, e.g. 
the rather large value of \g\ = 0.59 as recently determined. 

We have checked the sensitivity of our results to various parameters we used. The 
uncertainty in the strong coupling g = 0.59 ± 0.07, may change our results for direct 
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branching ratios by at most 15%. On the other hand, the uncertainty in A and A', and 
changing of sign of ai^2 is comparably neghgible. As to the vahics of //j, if wc vary it 
by a reasonable amount we can induce changes in the direct amplitudes by a few tens of 
percent. Althogeter, we feel that with the assumed calculational framework, the calculated 
direct amplitudes do not have uncertainties of more than 50%. Another feature which 
we neglect is the role of possible axial and scalar resonances, where the experimental 
situation is rather unsettled. 

As we explained in the text, the results (35) -(38) are obtained by using the experi- 
mental D — > Ktt amplitudes to calculate the inner bremsstrahlung. If we use the model 
for doing it, we get the result exhibited in (39)- (42), which do not differ from (35), (36), 
i.e. the — > K'^Ti'^'y decay, but are larger by a factor of about 2 in the amplitude in 
the case of K~7r~^'y decay. This is apparently related to the known difficulty of 

calculating the K~t[^ amplitude in the factorization approximation; therefore, we 

consider the results given in (37), (38) to be on a safer ground. 

If we disregard the contribution of vector mesons to the direct part of the radiative 
decays, the parity-conserving part of the amplitude is considerably decreased, by one 
order of magnitude in the rate in K'^Ti'^'-f decay and by two orders of magnitude 

m 7r+7. On the other hand, their contribution is not felt in a significant way in 

in the parity - violating part of the amplitudes. In any case, the detection of the direct 
part of these decays at the predicted rates, will constitute a proof of the important role 
of the hght vector mesons. The contribution of the vector mesons in the crossed channels 
is evident in the Dalitz plots of Fig. 7. We point out that the contribution of the vector 
mesons in Fi, as evidenced from the relevant graphs, is wholy determined by the A, A' 
couplings. 

Figs. 5a, 5b give the expected spectra for the direct component, which would be 
detectable in the region of high photon energioes. For the decay, the direct electric 
and magnetic transitions are of comparable strength. This prediction of the model should 
be testable, as it shiffts the peak of the spectrum to = 480 MeV, while if the magnetic 
transition is dominant it should peak — 400 MeV. In the — > K~7r'^'y decay 
the magnetic and electric components are likewise of nearly equal size, again testable 
in the spectrum. It is worthwhile to point out that the relative values of the parity - 
conserving amplitudes is rather large \A{D'^ K~n'^^)\'^ /\A(D~^ K^n^^)]^ ^ 7. The 
main reason for it are certain contributions (like -D?i), which appear in decay but 
are doubly Cabibbo forbidded in the decay. Also, as we pointed out, the A/ = 1/2 
amplitude is larger than the A/ = 3/2 one in the D Ktt channels. 

Finally, we wish to emphasize a most interesting implication of our calculation. When 
one compares the results obtained here for the radiative decays to nonresonant Ktt, with 
those previously obtained for the D — > K*^ [4,6,7,16], it emerges that the nonresonant 
channel is the more frequent one. Thus, while one expects for D'^ K*^'j a branching 
ratio (BR) of about 10~^ (this decay being doubly Cabibbo suppressed), the direct decay 
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— s> ^°7r+7 is expected to have a BR of ~ 3 x 10~^ in our model. To this one 
should add the IB component, which brings its BR to about 4 x 10~^ for /c > 50 MeV. 
The radiative decay K*^'~f is expected with BR of 0.5 x 10^'*. The nonresonant 

direct K~7r~^'~f decay we investigated here, has a BR of ^ 3 x 10~^ and including 

the IB component will occur with a rate 8 x 10~^ for k > 50 MeV. The experimental 
verification of this systematics will provide a check for the suitabilty of the theoretical 
methods employed. We should remark,however, at this point that we did not address the 
possibility of the decays D —>■ R + 'j , where R is a higher Kn or Ktttt resonance. To our 
knowledge, there is no calculation available on this topic. Our expectation is that such 
modes are at most comparable in strength to the D —>■ K*^ decay; preliminary data from 
BELLE indicate that this is the case in B decays. 

There is another feature of interest which we mention. Since there is interference 
between the direct electric and the bremsstrahlung component of the amplitude, these 
decays can also be used to test CP - violating effects in the amplitudes. 

We conclude by expressing the hope that the interesting features which these decays 
provide and were analyzed in this paper, will bring to an experimental search in the near 
future. 
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APPENDIX A: THE DECAY AMPLITUDES FOR D+ K'^n+^f 

Here we give the expressions for the sum of amplitudes in each row presented in Fig 1. 
The contributions which arise due to Oi operator are: 

.+ . fof^iy ■ q + v ■k)p-k 
A{ = —le- 



Jk V -k 



jI+ = _i / fpsfn p ■ q - {v ■ q){v ■ p) + V ■ k{M - v ■ p) p ■ k 
^ ~ M fK ^ v-p + A vk 

= ieg\ 



^sfDsfAvk){q-k){p-k) 
M fx V ■ k + vp + A 



■1 r2V-^A,i^-i%)U^f2A'-^A«.^^ 



V ■ k + A \ 2 ^ \ 3m2, ml 11 v ■ p + A \ 3 ^ m | 
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A2 = i\/2fDsf7r9K0'K''^9v>^eJMsM 



{v ■k){q- kf 



{v-{p + k) + A){ip + kf - + iVK^rriK*) ' 



The contributions coming from O2 operator are: 

r<+- ■jDfK{v-p){p-k), 



V ■ k 



^+ _ . fpfx p-q-{v ■p){v -q) + v ■ k{M -V ■p)p-k 



C+ = -ieg 



-1 



V ■ q + V ■ k + A 

fpfK {v-k){q-k) 
Jt^ P ■ k{v ■ k + V ■ q + A) 

qp\\ 1 \/2 



= -iV2fDfKgpn^gv>^eM- 



v-q + A 2 

{v ■ k){p ■ ky 



V ■ k 



+ 



3m^ m^^ 



V ■ {q + k) + A)((g + kf - m? + iTprUp) 



Next we give the expressions for the sum of amphtudes in each row presented in Fig 2. 
The contributions arising from the Oi operator: 



9- 



fosVM's V ■ q 



1 1 

+ 



Jk \v-k + A f^^/M V -p + v -kyv ■k + A vp + A^ 

ttiM — ■ q 



{p + kY - m|-^ - iVK*mK* 
2. 



-egp.,gp^VM ^''^^ 



'fx {k + qy - mp^ + iVprUp ' 

gu> _ qp\ 1 ' 

V^Vx'^""' \2>ml ml) (vk + A) 



1 

-T^-^Xgy 



MfD + x^MMJds 



V ■ q 



V • k + V ■ p 



I ^ /t\^i\4- fDsfiT . Q4> 



V ■ q 



fx "ml {v ■p + A){v k + v -p) ' 



The operator O2 gives the foUowing contributions: 



+ 9-, 



V • p 



\v-k + A {v q + v ■ k){v ■ k + A) J ' 
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1 + ^ 



, 1 ,^ JdIk . 
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V ■ k + V ■ q 
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APPENDIX B: THE DECAY AMPLITUDES FOR K'TT+^f 



The expressions for the sum of amphtudes (the Oi operator) in each row exhibited in Fig 
3. are: 

A^^^-iMe^{v-q + vk) , 

JK 



Al = tejMM, 



jDsfiT (p-k{v-p-M) p ■ k{p ■ q - {v ■ p){v ■ q)) 



fx \M{vp + A) M{vp + A){vp + v-k + A) 

Mq-k- Mp ■ q + M{v ■ q){v ■ p) - {v ■ q){q ■ k)\ 
M{v -p + v-k + A) y ' 

Z\ - ^Ag^—T 2 ^9v o„2 + ™2 



fosf-K {v ■ k){p ■ k){q ■ k) 



fx {v ■ k + v ■p + A) 



+ 



\ 



I 



A\ = i\/2fDsU9K*K^9v>^eJMsM 



{v-p + A) ' vk + A 

{vk){q-kY 
{v-{p + k) + A){{p + kf - m\. + iVK*mK*) 



The sum of amphtudes coming from operator O2, shown in Fig. 3 as C°, is vanishing. 

Next we present the expressions for the sums of amphtudes in each row shown in Fig. 
4. The results for the operator Oi are: 

jK {v-p + A){vk + v-p)\ 3 mi J 



-\/2 Jk {v-k-\-A)\ \3m5 ml) V M \3mJ m'lj (v-k + v-p) 
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aiM — ■ q 



{p + kf 



m 



K* 



fx {k + qY -ml + iVpiTLp 



Finally, the sums of amplitudes in each row due to the operator O^- 



a/2 ^ m?p)v-k + /S.\ {p + qY — ) ' 

r,^ = --f f ( dK-^K-r gpnj ^ 

' r''^'' gK^ {M^-m\,) \{p + kY-ml.+imK*TK*)^ {q + kY-m1 + impVp)^ 
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Figure 1: Feynman diagrams contributing to the formfactor Fi of the — > K^t:^^ 

decay. Diagrams denoted by Af^j {C^j) come from the operator Oi (O2). Sum of the 

contributions of each row is gauge invariant. In diagrams ^3^1, ^3^2; ^3,1 ^3,2 the 
photon couples to the heavy mesons with strength A'. 
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Figure 2: Feynman diagrams contributing to the formf actor F2 of the — > K^tt'^'j 
decay. Diagrams denoted by B^j (D^j) come from the operator Oi (O2). 
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Figure 3: Feynman diagrams contributing to the formfactor Fi of the D° — > ir~7r+7 decay. 
Diagrams denoted by A^^ (C^j) come from the operator Oi (O2). Sum of the contributions 
of each row is gauge invariant. In diagrams A^ ^ and 2 the photon couples to the heavy 
mesons with strength A'. 
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Figure 4: Feynman diagrams contributing to the formf actor F2 of the — > 7r"'"7 
decay. Diagrams denoted by B^j i^tj) come from the operator Oi (O2). 
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Figure 5: 



1 dr 



for the decay K°7r+7 (left) and ^ /C"7r+7 (right). Above: 



^total dm\2 

Direct parity - conserving (dashed hne) and parity - violating (putting Fq = 0) (full line) 
terms. Middle: tt-^ — -j^^ with T containing the full decay amplitudes, for model (dashed 

-L total df^l2 

line) and model+exp. (full line). For the latter, maximal Fq, Fi interference is exhibited. 

with r for the radiative decay calculated from model+exp. (full - line) 



1 dr 



Below: ^ , 

1^ total dmi2 

compared to pure bremsstrahlung emission (dashed - line). 



24 




0.£ 0.4 0.6 O.S 

Photon energy cut [GeV] 



[xlO-*] 

13r ^ 



12.5 



7.5 



£.5 




0.£ 0.4 0.6 O.S 

Photon energy cut [GeV] 



Figure 6: Branching ratio of radiative decays. Left: Decay — > K^n'^'y Right: De- 
cay — > ir~7r+7 Pull hne: Bremsstrahlung only. Long dashed line: all contributions 
included and pozitive Fq, Fi interference; short dashed line: negative interference 
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Figure 7: Dalitz plot of the parity conserving (above) and parity violating part (below) of 
the — > ir~7r+7 (left) and — > X°7r+7 decay (right). With gray levels on contour plot 
(left) and on z axis on the 3D plot (right) we present -^dV / {dmi2dm2z) in the logarithmic 

scale. Invariant mass mi2 = {P — ky is plotted on the x axis and 77223 — \J{P ~ pY on 
the y axis of contour plot. The x and y axes of the 3D plot are labeled. 
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